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PART  I 


Eq.  (2.9a)  reads  * 

Eq.  (2.9b)  reads  =  i(l  +  €j^/Ej^) 

^  ^  ' 

Second  line  of  Eq.  (3.I9I')  reads 
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Line  4  reads  "...  irrent  density  the  last  term" 


Eq.  (5*21)  reads 
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Eq.  (5*22)  reads 


R^(a5)  =  2^*^  ^^4  y'  ^ 


Reference  l4  is  A195»  556  (194-9  )• 
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Line  11  reads  "acting  within  the  energy  beuid..." 
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Eq.  (2.6b)  reads 


Eq.  (2.7)  reads 
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The  denominator  of  Eq.  (2.10)  reads 


.  2  . 

■0  „  2 


(e  +  X)2  +  a^  +iTi 
P  P  P 


The  denominator  of  the  integrand  of  Eq.  (2.15)  reads 


■n  * 


Line  8  reads  "Therefore,  the  i^l  term  may..." 


Reference  12  is  (The  Physical  Society,  London,  1959) 
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Part  1 

Excltons  and  PLasmons  In  Superconductors 
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I.  INTRODUtJTION 

In  the  original  theory  of  Bardeen,  Cooper,  and  Schrieffer^ 
an  approximation  to  the  ground-state  wave  function  of  a  supercon¬ 
ductor  was  obtained  by  a  variational  calculation.  Basic  to  the 

2/ 

theory  is  Cooper's  result—'  that  if  a  net  attraction  exists  between 
the  particles,  the  Fermi  sea  is  unstable  with  respect  to  the  forma¬ 
tion  of  bound  pairs.  The  BCS  ground-state  wave  function  is  formed 
from  a  linear  combination  of  normal  state-like  configurations  in 
whj-ch  particles  are  excited  to  states  of  low  energy  above  the  Fermi 
surface.  In  all  of  these  normal  configurations,  the  single -particle 
states  are  occupied  in  pairs  (I?f,  ),  so  that  interactions  other 

than  those  between  pairs  of  electrons  of  zero  net  momentum  and  spin 
are  neglected.  The  theory  leads  to  the  single  quasi -particle  exci- 
''dtion  spectrum  given  by  ■ 
energy  measured  with  respect  to  the  Fermi  level  emd  ^  is  the  energy 
gap;  that  is,  2^  represents  the  minimum  energy  required  to  excite 
a  pair  of  quasi-i>articles  from  the  ground  state.  The  quasi -particle 
excitations  are  fermions  and  no  boson  excitations  appear  other  than 
the  phonons. 

Ihis  independent  quasi-particle  approximation  has  been  sur¬ 
prisingly  successful  in  e]q>lainlng  the  thermodynamic  properties  as  well 
as  the  acoustic  and  electromagnetic  absorption,  the  nuclear  spin 
relaxation,  and  the  Meissner  effect  observed  in  the  superconducting 
state.  The  derivation  of  the  last  has  been  criticized  because  it 


(V 


where 


is  the  Bloch 
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is  not  strictly  gauge -invariant.  Biis  is  primarily  di«  to  the  neglec 
of  residual  Interactions  between  particles  in  states  and  ^  1?. 
These  interactions  give  rise  to  a  set  of  collective  excitations 
(hosons)  and  lead  to  a  gauge -invariant  description  of  the  Meissner 
effect. 

For  the  investigation  of  these  collective  excitations, 

.  3  / 

Anderson—'  euid  Bogoliubov,  Tolmachev,  emd  Shirkov—'  have  used  a 
generalized  time -dependent  self-consistent  field  or  random-phase 
approximation  (RPA).  Their  work  shows  that  in  the  superconducting 
shate,  the  plasmon  frequency  and  the  plasmon  coordinate  in  the  long- 
wavelength  limit  are  essentially  the  same  as  in  the  normal  state. 
They  have  also  suggested  the  existence  of  the  exciton  modes  lying 
within  the  energy  gap  which  we  investigate  in  the  main  body  of  this 
paper.  A  thorough  discussion  of  the  generalized  RPA  has  been  given 
hy  Rickayzen,-^  who  used  it  to  derive  the  complex  dielectric  constant 
of  a  superconductor  and  the  Meissner  effect  in  a  gauge -invariant 
manner.  The  BCS  quasi -particle  states  |a^  and  |  ^^  do  not 
satisfy  the  continuity  equation;  that  is,  |  V  *7  +  ^0. 

When  collective  modes  are  inclxxded,  the  current  suid  charge  density 
operators  7  and  p  are  decomposed  into  a  sum  of  individual-particle 
operators  and  collective  operators.  A  virtual  cloud  of  plasmons 
surrounds  each  quasi -particle,  producing  a  back- flow  current  which 
leeuls  to  over-all  charge  conservation  of  the  excitation.  Therefore, 
the  continuity  equation  is  satisfied  within  the  generalized  RPA. 
'nils  condition  is  sufficient  to  gueu-antee  a  gauge-invariant  form  of 
the  electromagnetic  response  kernel. 
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In  this  thesis  ve  interpret  the  exciton  mcxie  in  the  super¬ 
conductor  as  a  bovind  pair  of  quasi-particles  whose  center-of-mass 

propagates  with  momentum  The  exciton  spectrum 

is  investigated  through  the  generalized  RPA  equations  of  motion  pro¬ 
posed  by  Anderson  in  the  form  introduced  by  Rickayzen  involving  the 


v§/ 


,2/ 


quasi -particle  operators  of  Bogoliubov^  and  Valatin*^'  rather 

than  c^,  the  usual  electron  operators •  In  these  equations  we  make  an 
expansion  of  the  interaction  potential  in  terms  of  spherical 

harmonics.  It  is  found  that  excitons  may  be  characterized  by  the 
approximate  quantum  numbers  L  and  M  describing  the  symmetry  of  the 
states  with  respect  to  the  relative  coordinate  The  existence 

of  an  L-state  exciton  (corresponding  to  the  p,d,f,...  excitons)  is 
dependent  on  being  negative,  where  is  the  L-wave  part  of 
The  plasmon  state  corresponds  to  an  s-state  exciton  whose  energy  is 
greatly  increased  by  the  long-range  Coulomb  interaction. 

To  obtain  solutions  to  the  Anderson -Rickayzen  equations,  we 
take  matrix  elements  of  the  equations  between  a  state  with  one  col¬ 
lective  excitation  and  the  ground  state  which  has  been  renormalized 
so  as  to  include  the  zero-point  motion  of  the  collective  modes.  The 
results  give  two  sets  of  solutions  which  correspond 

to  what  Anderson  has  termed  odd  and  even  solutions.  We  show  that 
the  modes  are  lonphysical  and  that  the  modes  correspond 

to  the  exciton  states.  Ttie  quantum  numbers  L  and  M  are  found  to  be 
exact  in  the  limit  of  zero  center-of-mass  momentum^.  For  larger 
states  of  different  L  are  mixed,  although  the  mixing  is  small  for 
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«  1,  where  is  the  coherence  length.  The  magnetic  quantum 
number  M,  however,  remains  a  good  quantum  number  for  all  *5  if  the 
potential  has  no  crystalline  anisotropy.  The  exciton  energy  for  the 
^  =  0  case  is  plotted  as  a  function  of  the  L-wave  coupling  constant 
defined  by  g^  =  -  N(o)V^/47t,  where  N(o)  is  the  density  of  states 
in  the  normal  phase  at  the  Fermi  surface.  For  g^  >  g^,  the  excita¬ 
tion  energy  proves  to  be  imaginary  and  the  implications  of  this  with 
respect  to  the  original  BCS  ground  state  are  discussed.  The  M  3^  0 
excitons  may  be  considered  as  transverse  collective  excitations  since 
they  do  not  couple  with  a  longitudinal  field.  In  the  general  case, 
if  the  ground  state  is  formed  from  pairs,  the  exciton 

becomes  the  plasma  oscillation. 

In  Sec.  II  we  discuss  the  generalized  RPA  from  a  diagrammatic 
point  of  view.  Solutions  for  the  collective  excitations  are  obtained 
in  Sec.  III. 

Tsuneto^^  has  applied  Rickayzen's  cmalysis  to  the  problem  of 
the  surface  impedance.  While  he  finds  the  existence  of  a  precursor 
absorption  for  frequencies  within  the  gap,  his  results  give  an 
absorption  due  to  the  exciton  states  which  is  an  order  of  magnitude 

7/ 

less  than  that  observed  by  Ginsberg,  Richards  and  Tinkham^'  in  lead 


and  mercury. 
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In  Sec.  IV  we  consider  corrections  to  the  Anderson-Rickayzen 
equations  which  lead  to  a  new  type  of  exciton  of  a  particle-hole 
nature  closely  related  to  exciton  states  occurring  in  insulators. 

A  calculation  similar  to  Tsuneto’s  is  performed  for  the  electro¬ 
magnetic  absorption  due  to  these  new  exciton  states  in  Sec.  V,  As 
in  Tsuneto's  work,  the  ratio  of  the  surface  resistance  due  to 
excitons  to  that  of  normal  metals  in  the  extreme  anomalous  limit 
turns  out  to  be  about  an  order  of  magnitude  too  small  to  e^qplain 
the  observed  data. 


II.  EQUATIONS  OF  MOTION 

We  consider  a  system  of  electrons  interacting  via  an  effective 
two-body  potential  V,  whose  matrix  elements  in  the  Bloch  state  repre¬ 
sentation  are  given  by 


(H'X  ‘  VI  i  f  V- 


(2.1) 

This  potential  arises  from  both  Coulomb  and  phonon  interactions 
between  electrons  and  will  be  discussed  in  detail  below.  Ttie 
Hamiltonieui  is  e)q)res8ed  in  the  Heisenberg  representation  in  terms 


t 

of  the  operators 
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and  c,  „  which  create  and  annihilate  electrons 
in  Bloch  states  of  momentum  h  and  spin  a.  They  satisfy  the  usual 
Fermi  anticommutation  relations.  The  single -particle  Bloch  energies 
measured  relative  to  the  Fermi  energy  E^,  are  assumed  to  he  of 
the  form  (-ft  k  /2m)  -  E^,.  The  Hamiltonian  of  the  system  is  given  by 


(2.2) 


In  the  generalized  RPA  one  studies  the  time  evolution  of 
bilinear  operators  of  the  form 


C^)  -  ^  (2.3a) 

(2.5=) 


which  create  excitations  with  a  fixed  total  momentum  41$.  It  is 
helpful  to  consider  the  full-time  development  of  these  operators 
as  being  built  up  from  the  infinitesimal  change  of  the 
in  a  time  interval  6t;  for  exaiig)le^ 


operators 
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-t+ft)  -  i^(p  t) 

'  l>^(f^t)l  Si,  (: 

In  the  absence  of  the  interaction  V,  the  commutator  reduces  to 
^  ■» 

(?,t)  so  that  except  for  a  phase  factor,  the  operators 
are  independent  of  time.  We  call  any  operator  p/  an  eigenoperator 
if  its  time  dependence  is  given  simply  by  a  phase  factor.  The  equa¬ 
tion  of  motion 


L  HjyU>(  J  - (2.; 

for  the  operator  guarantees  that  when  applied  to  an  eigenstate 

I  of  H,  creates  an  eigenstate  |  of  H  with  an  excitation 

energy  From  the  Hermitian  conjugate  of  (2.5)  it  follows  that 

has  the  inverse  effect  of  p^.  That  is,  while  pj^”  adds  energy  to 
the  system,  subtracts  energy,  so  that  p^  and  p^  may  be  thoxight  of 
as  creation  aind  emnihilation  operators  of  excitations  of  the  system. 

A  knowledge  of  the  eigenoperators  and  their  eigenenergies  allows  one 
to  calculate  dynamic  properties  of  the  system  as  well  as  the  thermo¬ 
dynamic  functions. 

In  certain  cases  the  state  vanish  identically; 

for  example,  if  p^^"  creates  pairs  of  fermions  in  states  already 
occupied  in  I  p>  .  Another  example  is  if  the  operator  scatters 
excitations  already  present  in  the  initial  state,  in  which  case 
vanishes  when  applied  to  the  ground  state.  Both  cases  will  be  dealt 


with  in  the  next  section. 
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In  the  presence  of  the  interaction  V,  the  commutator  (2.4) 
is  complicated  by  the  presence  of  terms  involving  four  single- 
particle  operators  (c  and  c*^’s).  Therefore^  the  bilinear  operators 
b"^,  b,  and  p  are  no  longer  eigenoperators  of  H  and  one  must  include 
products  of  four,  six,  etc.,  single  operators  to  form  the 

in  this  case.  The  question  arises  whether  there  is  a  consistent 
approximation  in  which  the  eigenoperators  are  represented  as  linear 
combinations  of  the  bilinear  operators  b*^,  b,  and  P  alone.  Consider 
a  typical  term  in  the  commutator  arising  from  the  interaction  potential 

A'f' 

(2.6) 

This  expression  is  shown  in  diagramnatic  form  in  Fig.  1.  In  the 
diagram,  time  is  increasing  from  right  to  left  with  the  incoming 
particles  in  states”^  +  and  'Ir  entering  from  the  right.  The 
first  term  on  the  right-hand  side  of  (2.6)  is  represented  by  Pig.  l(a) 
in  which  the  interaction,  represented  by  a  dashed  line,  scatters  the 
spin-up  incoming  particle  to*^  +"5  +  q'  't  ,  creating  a  particle  and 
a  hole  in  states  -3'  'Ir  and  4r  ,  respectively.  In  Pig.  1(b) 


'f  i. 
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the  analogous  process  for  the  spir-down  particle  given  by  the  second 
term  in  (2.6)  is  shown.  If  at  time  t  -  o  a  pair  of  single  particles 
is  excited,  at  time  5t  there  is  a  finite  probability  that  a  particle- 
hole  pair  has  been  created  from  the  background  of  particles  in  the 
Fermi  sea,  with  the  incoming  particles  scattering  to  new  states.  In 
the  next  interval  of  time  a  similar  process  may  occur  involving  any 
of  the  four  excitations,  and  in  general  the  **bare"  incoming  particles 
will  create  a  complicated  cascade  of  excitations  leading  to  a  decay 
of  the  initial  state.  In  the  generalized  random  phase  approximation 
one  keeps  only  those  terms  in  the  commutator  which  conserve  the  number 
of  excitations  allowing  for  both  forward  and  backward  propagation 
in  time  (see  below).  This  procedure  corresponds  to  a  linearization 
of  the  equations  of  motion  by  replacing  two  single-particle  operators 
in  each  term  by  a  c-number  given  by  the  expectation  value  of  this 
pair  of  operators  with  respect  to  a  fixed  state.  If  this  state  is 
chosen  to  be  the  BCS  ground  state,  defined  by 

\p-J\LUi,tV^bl(pnQ)^ 

where  |  0^  is  the  state  with  no  particles  present,  conservation 
of  momentum  and  spin  leads  to  nonzero  average  values  only  for  the  op¬ 
erators  b^  (o),  ^j^(o),  and  S  In  terms  of  the  pareuneters 

u^  and  Vj^,  these  averages  are 
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(2.8b) 


The  parameters  and  are  given  by 

-h  ( \  )"**  j 

z^  =  +  (f/-4/£R)'^, 


vhere 


and  ^  satisfies 


(2.9a) 


(2.9b) 


(2.10) 


(2.11) 


This  prescription  gives  a  unique  linearization  of  the  equations  of 
motion  since  for  q  ^  0  there  is  at  most  one  pair  of  operators  with 
zero  total  momentum  and  spin  in  each  term.  Ihe  terms  retained  within 
this  approximation  are  shown  in  Fig.  2. 

(1)  As  shown  in  Pig.  2(a),  the  conventional  particle -particle 
scattering  vertex  arises  from  the  first  term  in  (2.6)  when^‘  * 
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The  factor  of  ~  in  front  of  V  is  cancelled  by  the  term  in  the  inter¬ 
action  with  spins  opposite  to  those  in  (2.6).  This  cancellation  of 
the  factor  of  ^  occurs  in  each  vertex. 

(2)  Another  possibility,  shown  in  Figs.  2(b)  and  2(c),  is  for 

the  scattered  incoming  particle  to  enter  a  bound  state  with  the  other 

incoming  particle,  the  outgoing  excitations  being  the  particle-hole 

pair  created  from  the  sea.  This  possibility  is  allowed  for  in  the 

linearization  by  including  the  finite  average  ^ 

which  may  be  regarded  as  the  amplitude  for  the  pair  to  enter  the  q  «  0 

bound  state,  which  is  raacroscopically  occupied  in  l'^0>  Since 

a  finite  fraction  of  all  the  electrons  occupy  this  bound  state  in  the 

superconducting  state  (corresponding  to  the  finite  fraction  of 

helium  atoms  occupying  the  k  =  0  state  is  superfluid  He  ),  the  small 

1  /2 

fluctuation  ^  N  ^  in  the  number  of  pairs  N  described  by  (2,7) 
leads  to  no  difficulties  in  a  large  system.  Notice  that  in  Figs.  2(b) 
and  2(c),  the  incoming  pair  of  i>articles  is  transformed  into  a  particle- 
hole  pair  by  the  interaction.  Therefore,  b (^)  and  P,  (^  are  coupled 
in  the  equations  of  motion. 

(5)  In  addition,  there  is  the  possibility  that  the  scattered 
incoming  particle  enters  the  bound  state  with  the  pcurticle  created 
from  the  sea,  leaving  the  hoic  and  the  other  incoming  particles  as 
the  outgoing  excitations,  as  shown  in  Figs.  2(d)  and  2(e).  Due  to 
the  presence  of  the  bound  state,  the  incoming  spin-up  particle  in 
Fig.  2(d)  is  transformed  into  a  hole  in  the  state  of  opposite 
momentum  and  spin.  In  the  next  instant  of  time  the  inverse  process 


15 


may  occur.  It  is  clear  that  the  equations  of  motion  are  simplified 
if  one  introduces  “quasi -particle”  oi)erators  7^^  which  are  the 
proper  linear  combinations  of  particle  and  hole  creation  operators 
to  account  for  these  processes.  The  appropriate  transformation, 
introduced  first  by  Bogliubov  and  by  Valatin,  is 


^  ^  J  (2.12a) 

(2.12b) 


For  mathematical  simplicity  we  will  follow  Rickayzen  by  expressing 
the  final  linearized  equations  in  terms  of  quasi -particle  variables. 

(4)  The  exchange  contributions  to  the  single-particle  lines 
are  shown  in  Figs.  2(f)  and  2(g).  As  is  well  known,  they  lead  to  an 
anomalously  low  density  of  states  at  the  Fermi  stirface  in  the  normal 
metal  unless  a  screened  interaction  is  introduced.  This  point  is 
discussed  below.  The  exchange  self-energy  vertex  can  be  accounted 
for,  along  with  process  (5),  by  the  quasi-particle  transformation 

(2.12). 

(5)  Finally,  the  xonscattered  incoming  particle  may  enter  the 
bound  state  with  x^he  i>article  created  from  .ne  sea,  leaving  the  hole 
and  scattered  peurticle  as  the  outgoing  excitations,  as  shown  in 

Figs.  2(h)  and  2(i).  As  in  process  (2),  the  pair  of  incoming  particles 
is  tramsformed  into  a  pcurticle-hole  pair  by  the  interaction.  In  the 
limit  q  0,  process  (2)  is  more  important  than  (^)  in  forming  the 
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plasmon  state.  Since  the  momentum  transfer  is  always  in  the  former 

2  2 

process,  the  large  matrix  element  of  the  Coulomb  interaction  kite^/q 
dominates  the  latter  vertex  in  which  the  momentum  transfer  may 
assume  any  value.  Anderson  and  Rickayzen  have  neglected  processes 
(4)  and  (5),  suggesting  that  their  effect  is  primarily  to  renormalize 
the  single -particle  energies  and  the  effective  interaction. 

The  terms  occurring  in  the  linearized  equation  of  motion  for 

are  shown  in  Fig.  5  and  bear  a  close  resemblance  to  those 
shown  in  Fig.  2.  In  the  conventional  RPA  for  the  excitations  in  the 
normal  state,  only  the  polarization  vertex  C  Fig.  5(b)*]  is  retained. 
Tne  so-called  exchange  scattering  correction  shown  in  Fig.  5(a),  when 
ccanbined  with  the  polarization  vertex,  approximates  the  time  evolution 
of  by  graphs  of  the  type  shown  in  Fig.  4.  In  the  limit  q  0, 

the  exchange  correction  to  the  plasmon  frequency  vanishes.  Since 
matrix  elements  of  the  equations  of  motion  are  taken  with  respect  to 
RPA  eigenstates,  two  pairs  may  be  spontemeously  created  from  the 
vacuum  and  may  interact  with  the  incoming  excitations  as  in  Fig.  k. 
This  process  may  be  viewed  as  a  propagation  of  the  excitations 
backward  in  time,  familiar  in  the  Green *8  function  formulation  of  the 
problem. 

In  the  generalized  RPA  for  the  superconducting  state  the 
presence  of  the  bound  state  gives  rise  to  the  vertices  represented 
in  Figs.  5(c),  (d),  (g),  and  (h),  so  that  eui  incoming  particle-hole 
pair  can  be  transformed  into  either  a  pair  of  particles  or  a  pair  of 


15 


holes.  Therefore,  the  operators  (^)  and  coupled  by 

the  density  operator  vertices  occurring  in  the  time 

development  of  b  , C^)  are  identical  to  those  in  Pig.  2  except  that 

JX 

all  arrows  are  reversed  and  the  momentum  is  replaced  by 

We  turn  now  to  the  question  of  screening*  Within  the  random- 
phase  approximation  to  the  normal  state,  the  screened  interaction 
line  is  represented  in  the  limit  of  small  wave-vector  ^  by  a  sum 
of  diagrams  of  the  form  shown  in  Fig.  5*  Rickayzen  has  shown  that 
the  dielectric  constant  is  essentially  unaffected  by  the  pairing 
correlations  occurring  in  the  superconducting  state.  It  is  easily 
seen  that  the  vertices  2(b),  2(c),  and  5(^)  are  automatically  screened 
within  the  RPA  through  the  presence  of  the  polarization  vertex 
C  Fig.  5(b)  3  in  the  linearized  equations.  For  example,  when  the 
vertex  2(b)  is  followed  in  time  by  a  series  of  vertices  5(b),  the 
effect  is  to  replace  the  bare  interaction  line  in  2(b)  by  the 
screened  line  shown  in  Fig.  5*  Therefore  in  vertices  2(b),  2(c),  and 
5(b),  the  unscreened  interaction  must  be  used.  The  potential 
is  given  by 


(2.15) 


where  -K/i  is  the  energy  of  the  excitation  involved.  Also,  v  ^  is 

<1 

the  bare  electron -phonon  interaction  matrix  element  introduced  by 
Bardeen  and  Pinea—^  and  is  the  bare  phonon  frequency.  It  is 


l6 


essential,  however,  to  introduce  the  interaction  screened  by  the 
dynamical  dielectric  constant  in  the  remaining  vertices  since  it  is 
impossible  to  replace  the  bare  interaction  line  by  the  screened  line 
through  an  iteration  of  vertices  occurring  in  the  linearized  equations. 
The  screened  potential  is  of  the  form 


J 


where  the  dynamical  dielectric  constant  is  given  by 


(2.14b) 


Here,  ^  ~  is  the  electronic  screening  wave 

number.  In  a  more  complete  treatment  involving  coupled  equations  of 

motion  for  the  electrons  euid  the  lattice,  the  energy  ^oi  would 

K,q 

presumably  be  given  in  terms  of  the  quasi-particle  excitation  energies. 

For  siii5)licity,  we  neglect  the  vertices  shown  in  Figs.  2(h), 
2(i),  and  3(a)-  We  also  neglect  the  exchange  self -energy  correction 
since  it  simply  renormalizes  the  single -particle  energies.  With 
these  approximations,  one  obtains  the  equations  first  given  by 
Anderson : 


(2.15a) 
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(2.15b) 


C^/^/^tf]  =  li  ff  (p^ 

I"  p 

''  '^Htf  21^  ^  (-'J' 'j  ’  ■’ 

In, p.^^4' (f)}  --(^-ii^)p^-^*(f)-(vC-^^)\l (pf  f 

"  A), i))^ (p +^teLj  Kp'^^+^22k^^  bt,' Cp 

-  TI-KVklL  l>K<HI  (-p  ■ 


(2.15d) 


The  density  operator  p(q)  is  given  by 

V- 


As  mentioned  above,  the  equations  can  be  considerably  simplified 
by  tremsforming  to  quasi-particle  variables.  The  Anderson-Rickayzen 
equations  are  then: 


18 


^  (ftfpf) 

(2.16a) 


(2.16c) 


The  coherence  factors  are  defined  by 

(2.17a) 

(2.17b) 

(2.17c) 

~  ^ (2 . 17d ) 


and  the  three  collective  variables  are 
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l^k  (p  =  L  l>k'  ~  l>k 


di^(p  =Z  W'-r, W  [i#<p  A'tjr  ('p-i 


(2.l8c) 


From  (2.16c)  we  see  that  half  of  the  normal  mode  operators  are 
of  the  form  ^'k.a*  eigenvalue  'Riese 

operators  describe  scattering  of  excitations  already  present  in  the 
Initial  state  and  vanish  when  applied  to  the  ground  state.  Since  we 
will  always  tcOce  matrix  elements  of  the  equations  of  motion  between 
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the  ground  state  and  an  excited  state,  these  quasi -particle  conserving 
operators  may  be  safely  neglected. 

III.  SOLUTIONS  OP  EQUATIONS  OP  MOTION 
Por  the  analysis  of  the  plasmon  and  exciton  modes  at  temperature 
T  *  0  we  begin  with  the  Anderson-Rickayzen  equations  of  motion  (2.l6) 
for  the  pair  operators  ^k+ql^kO'  “^^st  be  kept  in 

mind  that  the  equations  have  been  linearized  with  respect  to  the 
ground  state  involving  s-state  pairing  between  electrons  of  opposite 
spin  and  momentum,  as  our  results  depend  critically  upon  this  fact. 

The  collective  variables  defined  by  (2.18)  are  substituted  into  the 
equations  in  order  to  obtain  them  in  a  form  involving  only  the 
Bogoliubov-Valatin  quasi -particle  operators: 


(3.1a) 
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3£;>f.p2  yoi,mc^'.f)(fkl^]f^,  -/*-^,^-o)- 

'th  (^jf)Z-  +/k'+^<  i^'o). 

(5.  It) 

Those  operators  p.^  (^ )  are  now  considered  which  are  linear  combina¬ 
tions  of  the  bilinear  products  of  ^j^'s  appearing  in  the 

two  equations  of  motion  (5*l)>  and  which  create  one  elementary 
excitation  of  type  a.  Ihus  we  desire 

(5.2) 

with 


where 


lo> 


is  not  the  original  ground  state  of  BCS,  but  the 
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renormalized  ground  state  with  1  0^  =0*  Th®  quantity 

represents  the  energy  of  the  excitation  created  by  the  operator 
The  elementary  excitation  may  'be  any  one  of  the  three  types 

involved  in  the  theory:  a  pair  of  excited  quasi-peurticles  in  scattering 
states,  a  plasmon,  or  an  exciton. 

From  Eq.  (5*3)  and  the  discussion  of  Sec.  II,  we  have 

[H.wjct)]!  o>  =^£:^c?)K„'t'(t)io>  ■ 

Since  the  commutator  c  is  related  to  the  time  derivative 

of  (^),  the  matrix  element  of  (^)  between  the  ground  state  (  0^ 
and  the  state  |  l(^,a)^  containing  one  excitation  of  energy 
must  have  the  time  dependence  exp  [iXl^(t)t]  .  Now,  Eq, 

(5*2),  expresses  ,  within  the  RPA,  as  a  linear  combination  of 

the  bilinear  products  and  7^+q]^7j^0^  that  we  may  write 

the  inverse  transformations  as 

H/-0  -- 1  p 

(5.1ft.) 

■DEiking  matrix  elements  of  Eq.  (5.*+)  between  |  0  ^  euid  I  l(9,a)> 


and  using  the  orthonormality  property  of  the  excited  states,  we  find 
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(3.5a) 

The  solution  for  the  exciton  mode  dispersion  relation  is 
dependent  on  taking  matrix  elements  of  the  equations  of  motion  (3-l) 
between  the  states  1  o>  and  1  l(^,a)  ^  and  using  the  relations 
(5,5a.)  and  (5*5b)  so  that  we  obtain  a  set  of  c-number  equations, 

The  resultant  system  of  linear  equations  may  then  be  solved  for  the 
normal  mode  frequencies  and  the  transformation  coefficients  f  and  g. 

By  taking  matrix  elements  of  (5«l)  we  obtain ^ 

K, 

ihCK.p  L  KiP,  W  no^.p  ifc 


(3.6a) 


2k 


A 

iw.p  Z  V(lPM(^'^pljk(l?.f)-^(l?'pl 

(5.6b) 

From  (5*6)  it  is  evident  that  an  explicit  form  for  VCS^^*)  must  be 
chosen  in  order  to  proceed  further c  As  emphasized  in  the  foregoing, 
the  BCS  ground  state  about  which  the  Anderson-Rickayzen  equations 
have  been  linearized  is  one  involving  s-state  pairing.  Thus  in  the 
absence  of  crystalline  anisotropy,  the  q  -*  0  8olu":ions  must  transform 
according  to  the  irreducible  representations  of  the  full  rotation 
group,  i.e.,  the  spherical  harmonicsc  Because  of  this  fact,  we 
expand  the  two  body  potential  V(?,^*)  in  terms  of  spherical  harmonics. 
Ttie  coordinate  system  is  chosen  so  that^  lies  along  the  polar  axis 
with  0  and  cp  the  polar  and  azimuthal  angles  of  the  wave  vector 
and  6  and  ®  the  analogous  quantities  for  If  ^  is  the  angle  between 

k*  and  k,  the  use  of  the  addition  theorem  gives 

0* 

00  '' 
xso  M0-A.  ■'  ' 


(3.7) 
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where 

V^(k,k')  =  (W2^  +  1)^  V^(k,k')  . 

A  further  approximation  is  made  in  setting  V^(k,k')  =  a  nonzero 

constant,  for  I  €]  and  zero  otherwise.  The  quantity^ 

^  *  c  ^ 

is  the  average  phonon  energy  of  the  order  of  the  Debye  energy.  We 
define  the  coupling  constant  gj^  by 

=  -N(0)  Vj^Ait  .  (5.8) 

The  BCS  coupling  constant  is  related  to  g^  by 

EO  =  -N(0)  vyi*,  .  N(0)  Vjj,3  >  0  . 

It  is  convenient  to  introduce  three  new  variables 

r^(f)  --ZuVfWj. 

where  the  subscript  a  has  been  dropped  from  both  sides  of  the  equa¬ 
tions  for  simplicity.  Equations  (3.6)  then  express  the  transformation 
coefficients  f  ewid  g  in  terms  of  the  new  variables  A,  P,  and  Z.  By 
substituting  these  expression  into  the  defining  relationships  (5*9)> 
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we  obtain  the  following  coupled  integral  equations  to  determine  the 
e i ge  nf re  que  nc  y  a  (q): 


(Ij)  [^j2 (pf_  jjK 


(3.10b) 
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\  ziipM<^,p1(p* 

<^SL(pMlp7-  h»<  ^*' 


tlons  It  1=  l^^dl^tely  se,t  ttot  one  good  qnnntnm 
'^7:  nn  e^ttetton  ts  tne  ..ette  ,^ntn„ 

;  7  ne^ntnee  «  -  «. 

>,o-r  M  In  the  sum  over  R,  the  angui^ 

'  .  uantities  involved  are  the  spherical  har¬ 
as  the  only  (p-dependent  quantit  center-of- 

.onics.  Thus,  M  is  a  good  quantum  number  regardless  o 

Qiass  momentum  Aq* 

(1)  q  '  0 

e.„ter-of-mans  mo»entn«,  E,s.  D-IO)  81™ 
In  the  case  of  zero  center  of  ma 

«>,or-  This  follows  since  neither 

L  as  an  additional  good  quantum  n  • 

the  energy  quasi -particle  pa 

the  coherence  factors  nor  angular  part  of 

dependent  on  tn.  podnn  an.te  U  tnt.  case.  Tne  an. 

the  sum  \  then  reduces  to 


The  sum  is  converted  into  an  integral  by  letting 

where  the  volume  v  of  the  normalization  "box  is  taken  as  unity.  The 
radial  integrals  over  k  are  all  of  the  form 


where  each  of  the  quantities  a,b,c,  ...  is  one  of  the  coherence 
factors,  the  energy  (o)  of  the  independent  quasi-particles,  or  the 

A 

excitation  energy /ft/1.  The  integration  over  the  magnitude  of  Vis 
replaced  by  an  integration  over  the  Bloch  state  energy  as  meas¬ 
ured  from  the  Fermi  surface,  by  setting 

N(o)<li-j  (5.12) 

where  we  have  made  the  approximation  of  a  constsuit  density  of  states. 
The  approximation  leads  to  an  error  of  order  Ajo^/Ej,  =  lo”^.  The 
Integrals  0  are  only  performed  over  the  region  <  e 

since  the  potentials  have  been  set  equal  to  zero  outside  this 
energy  band.  Using  (5.12),  Eqs.  (5.10a)  and  (5.10b)  for  the  q  -*  0 


case  are  written  as 
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(l  -  \4 

^0, 

=  / 

f 

From  these  equations  it  is  seen  that  the  direct  Coulomb  interaction 
Une  /q  involved  in  Z(^)  only  appears  for  the  L  =  M  =  0  state.  It 
will  be  shown  below  that  this  state  has  a  solution  corresponding  to 
a  plasma  oscillation  with  the  usual  plasmon  energy 

'hiZp-'ft  loe^ 

and  lies  far  above  the  gap  2A  10  ^  ev.  In  this  section  only  the 
M  ^  0  cases  will  be  considered,  in  which  the  right-hand  sides  of 
Eqs.  (2.15)  become  zero.  Since  the  integrand  of  is  odd 

about  the  Fermi  surface  within  the  constant  density  of  states  approxi¬ 
mation,  vanishes  and  there  is  no  coupling  between  the  A  and 

r  modes.  The  excitation  energies  for  the  L  0  modes  with  zero 
center-of-mass  momentum  are  then  determined  by  the  conditions : 


(5.13b) 


50 


(i-v,r^j=o 

0  -  Vi  =  0, 


(Ajj^  mode)  ,  (5.1^) 

(r^mode)  .  (3.l4b) 

and  I  using 

Eqs.  (5- 13)  become: 


Setting  X  =  (‘ftX2./2A)  <  1  in  the  integrals  Ivn2^ 
the  definition  (5-8)  of  the  coupling  constant  g. 


( 


I  -  J. 


(5.15a) 


(5.15b) 

Values  of  X  «  (^II/2A)  are  plotted  as  a  function  of  the  left-hand 
sides  of  these  equations  in  Fig.  6.  The  plot  shows  that  when  g^^ 
becomes  larger  than  g^,  the  frequency  n  of  the  Pj^  mode  becomes 
imaginary,  indicating  that  the  system  is  unstable  when  described 
by  a  ground  state  formed  with  s-state  pairing.  Therefore,  if  g^^  is 
the  largest  coupling  constant  present,  the  groxind  state  should  be 
formed  from  pair  functions  having  L-type  symmetry.  The  pair  spin 
function  is  singlet  or  triplet  depending  on  whether  L  is  even  or 
odd,  since  the  wave  function  describing  the  exciton  state  must  be 
antisymmetric  on  the  interchange  of  all  coordinates  of  the  quasi- 
particle  pair  involved. 
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The  growth  of  the  P  modes  for  g  >  g  also  indicates  that 
LM  L  0 

the  modes  have  no  physical  existence •  As  is  seen  in  Fig.  6,  a 
cannot  exist  unless  >  g^.  However,  when  such  a 
coupling  strength  is  reached,  the  corresponding  exciton  is  un- 
Stable  so  that  the  system  decays  before  the  A  mode  can  come  into 
existence.  Figure  6  also  indicates  the  2L-fold  M  degeneracy  of 
q  =  0  L- state  excitons. 

It  should  be  mentioned  that  a  continuum  of  scattering  state 

solutions  is  obtained  from  (5*l^b)  corresponding  to  the  vanishing 

of  the  denominator  of  the  integrand.  One  such  state  exists  between 

two  successive  unperturbed  levels,  E  +  E  .  Although  the  energy 

K  K+q 

of  a  scattering  state  solution  is  unaltered  from  its  value  in  the 
absence  of  interactions,  its  wave  function  is  strongly  modified  since 
each  particle  is  surrounded  by  a  depletion  of  the  same  type  of 
particle  leading  to  the  backflow  picture  mentioned  above. 

(2)  ^  Finite  Case 

From  Eq.  (j.io)  it  is  seen  that  L  is  not  strictly  a  good 
quantum  number  for  the  case  of  finite since  the  coherence  factors 
and  now  have  a  polar  angle  dependence.  Because  of  the  com¬ 

plexity  of  this  dependence,  the  sum  cannot  be  carried  out  exactly. 
We  approximate 


(3.15c) 
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where  p  \x  =  cos0^  and  is  the  velocity  of  a  particle  at 

the  Fem^i  surface.  This  leads  to  an  error  of  order  «  1.  The 

integrals  are  of  the  same  form  as  those  in  the  q  =  0  case •  To 

perform  the  angular  integral^  we  expand  the  denominator  of  the 
integrand 


(5.16) 


in  powers  of  p.  This  procedure  is  valid  so  long  as  P  -  2A. 

The  integrals  over  k  are  then  of  the  form 


ab. 


=  I 


ab. . 


+  ^  ^ 
^  ab . . .  ab . 


(3^17) 


with  superscripts  indicating  the  powers  of  p  involved.  Keeping 

2 

terms  through  order  p  and  using  the  relations 

COS0  =  n  =  (W5)^^^  Y^q(0) 

and 

cos^e  =  =  |(W5)^^^  Y2q(0)  +  ( Yqo  , 

the  eqxiations  for  A  and  P  (5-10)  become 


Aui  (p  Xc  + 

J-*cuit  X  ^  )i/y  l/«  p)  i, 

(3.18a) 
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+ifhu^)Yoo  +f(f4i:MZ  Yy^Yj^r,.f]. 

(3.i8b) 


With  the  relation 


W(^AH)  J 


A 


where  the  C's  are  usual  Clehsch-Gordan  coefficients,—'*^  Eqs.  (j.lS) 
become 


* 
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r  laT 

L  iW(^LfO\ 


X  CCozLjoMM)  CiOiL  'jOOo)  AxhLP'^ 

JL 

X  (uljocoI/Uh  (p  +  Z  [1^^;]^ 

;<  (ijll;oMH)COjIIi  Doo)  >2^  (pi. 


(5.19a) 
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y H (LCoaL'  ooo)  Ctn  (p 


,  n  /V  /^\'V'  5~ r^iMtiili- 

t  (M  Lj  ooo)  Ijrf  (p  ■*-  -i-hajh  ^  L  ¥/rr^HjJ 

X  ^  UA  tj  oM^^)  CO  A  o  i>oo)Ajin  (pi. 


(3.19t) 
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As  in  the  q  0  case,  the  Coulomb  field  represented  by  the  presence 
of  the  Z(q)  term  does  not  couple  into  the  equations  of  motion  except 
for  the  longitudinal  modes  M  =  0.  Discussion  of  this  case  is  deferred 
and  the  transverse  cases  M  ^  0  are  now  considered.  For  a  given 
M  ^  0,  Eqs.  (3.19)  represent  a  set  of  2N  linear  simultaneous  equa¬ 
tions  in  A-.,  and  where  N  is  the  number  of  terms  present  in  the 

spherical  harmonic  decomposition  of  the  two-body  interaction  (3«T)« 

It  follows  that  for  a  given  set  of  V- *s  the  normal  mode  frequencies 
of  the  system  may  be  obtained  by  setting  the  determinant  of  the  coef¬ 
ficients  of  the  and  equal  to  zero.  Once  the  frequencies 

have  been  obtained,  the  A  *s,  F  *s,  and  the  transformation  coeffi- 
cients  f  and  g  may  be  determined. 

For  simplicity  we  consider  the  case  for  which  all  but  two 
of  the  vanish.  It  is  assumed  that  the  two-body  potential  con¬ 

sists  of  a  term  V^,  corresponding  to  the  BCS  parameter  and  another, 

Vj^,  representing  the  angular  dependence  of  the  interaction.  Since  M 
has  been  taken  as  nonzero,  it  is  seen  that  the  simplified  V  and  V 

U  JLi 

potential  allows  the  modes  to  be  characterized  by  a  quantum  number 

L  within  the  approximations  of  the  calculation,  due  to  A  and  F 

UN  CM 

vanishing  identically  for  M  ^  0.  Thus,  we  may  speak  of  a  p-,d-,... 
state  exciton  when  the  additional  term  in  the  potential  has  L  -  1,2,... 
type  angular  dependence. 

If  the  potential  contains  s-  and  p-wave  potentials, 

V(z  m-Vp  Y^(6.f)  Yco  (&.  + 


(3.20) 
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the  dispersion  relations  obtained  from  (5.19)  are  found  to  be 


^  (5.21a) 


X  =.  Cr^a.  f  M  T^a.)^  r  ±  1  (JJ  mo</e^] . 


(5.21b) 


We  discard  the  A  mode  since  it  does  not  exist  if  the  system  is 

stable.  The  dispersion  relation  (3*21b)  for  the  T  mode,  when 

1^+1 

rewritten  in  terms  of  e3q>licit  expressions  for  the  integrals 
and  becomes 


JL- 


£2. _ ff-  X*- _ 


(5.22) 


Where  x  <  1.  This  dispersion  relation  is  plotted  in 

Pig.  7  for  two  values  of  with  g^  »  0.25.  From  the  figure,  it  is 

seen  that  the  curve  intersects  the  origin  for  g^  =  g^.  For  a  valxie 

®1  ®0  minimum  value  of  x  =  given  by  (l/g^^  -  ^/Bq)  = 

2  1  /2 

X  arcsin  x  /(l  -  x  )  ,  in  aigreement  with  the  results  of  the  last 

m  m  m 

section  for  the  q  0  case. 
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(5)  The  s-State  Exciton 

The  above  discussion  was  restricted  to  that  of  the  trans¬ 
verse,  M  0,  excitations  in  which  the  Coulomb  interaction  term 
Z(t)  did  not  enter  into  the  equations  for  )  and  Be- 

fore  discussing  the  M  =  0  cases,  it  should  be  emphasized  that  the 
equations  of  motion  (5*1 )  which  are  the  basis  of  this  paper  are 
those  linearized  by  Anderson  about  the  BCS  ground  state  based  on 

3  12/ 

s-state  pairing  of  the  electrons*  As  Anderson^^-^ — '  has  pointed 
out,  it  is  the  s-state  exciton  which  corresponds  to  a  plasmon 
excitation,  due  to  Z  it)  coupling  into  the  equations  of  motion. 

The  L  =  0  mode  is  considered  in  the  q  0  limit.  Because 
of  the  singular  nature  of  the  direct  interaction,  it  is  not  possi¬ 
ble  to  set  q  S  0  in  the  calculation,  so  that  the  limit  q  0 
must  be  taken.  For  our  starting  point,  we  consider  Eq.  (3- 15b)  for 
the  ^qq(^)  mode  in  the  q  0  case: 

0  (3.25) 

From  the  definitions  (3-5)  and  (3-9)  an  expression  for  Z(q)  is  obtained 

Zc?) - - r  y 

(3.24) 
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Since  the  L  =  0  mode  excitation  energy  is  being  considered^  only 
the  term  in  (5.24)  need  be  used  in  substituting  for  Z(q) 

into  (5*25).  Rearrangement  of  terms  then  gives: 


'j-  lini 

Q^o  0 


Cr,;*+i'«-/3Xi-Kor«j  1. 


(3.25) 


p 

Since  V^(q.)  ~  l/<l  >  Eq.  (5.26)  indicates  that  in  order  for  the 
limit  to  be  finite,  the  terms  in  the  numerator  which  are  independent 
of  q  must  vanish: 


(3.26) 


The  validity  of  (3.26)  is  shown  by  considering  the  explicit  form 
of  the  integrals  involved: 


4o 


(5.27a) 


-  '^•/2 


•^0), 


m)A/£ 


dL 


'^0>, 


'"  v(t)  <■ 
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I  (MQ-K 

A(—L _ L 

£  aa)yi 


;  1::^  .  _ 

hsvi 


-  £A  f_J. 

*>-n.  ( 


(5.27b) 


where  the  BCS  integral  equation—^  for  has  been  used  to  obtain 
the  first  term  on  the  right.  With  the  use  of  these  relations, 
(5.26)  becomes : 


4l 


yc 


-h 


V  ^0. 

With  the  validity  of  (5-26)  established,  (5*25)  reduces  to 


/ 


im  _ 


(3.28) 

To  determine  the  existence  of  a  plasma  oscillation  for  the  L  «  0 

mode,  (5»28)  must  have  a  solution  for  x  «  ^jQ./2A)  »  1.  Under  this 

0 

condition  the  term  ^q^v/2  denominator  is  much  less  than 

unity  and  may  be  dropped.  'Rie  integrals  involved  in  (5.28)  are 
evaluated  for  x  »  1  so  that,  to  order  l/x  ,  (5-28)  reduces  to 

(5.29) 

Using  Vjj(q)  »  kne^ jc^  and  e^^^N(o)  »  (5/2jr^)(ha)p/2A)^,  where 
0)^^  ■  4ime^/m,  (5.29)  gives  JCI  »  cd^  so  that  the  excitation  fre¬ 
quency  of  this  mode  is  the  plasma  frequency. 
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(4)  The  L  =  1,  M  =  0  Mode 

To  complete  the  investigation  of  the  collective  states  present 
when  only  the  and  terms  are  kept  in  the  potential  expansion 
(5.7),  we  must  determine  the  dispersion  relation  for  the  r^Q(q)  mode. 
Setting  M  »  0  in  (5*19^)  ve  obtain  two  simultaneous  equations  involving 
roo(q)  and  There  is  no  mixing  of  these  modes  in  the  equations. 

The  "oo  dispersion  relation  gives  the  plasma  frequency  as  discussed 
above  while  the  mode  dispersion  relation  becomes 


(r;,.  ^  . it;,.). 


(3.50) 


In  Sec.  Ill  (2)  we  found  the  dispersion  relation  for  the  modes 

to  be 


(5.21b) 

Thus  the  dispersion  relation  can  be  obtained  by  letting 

<1  J  3  in  (5.22),  indicating  that  for  a  given  wave  vector"^  the 
excitation  energy  of  the  longitudinal  mode  is  raised  above  that 

of  the  transverse  modes. 

IV.  CC«RECTIONS  TO  THE  ANDERSON-RICKAYZEN  EQUATIONS 
We  consider  here  the  terms  in  the  linearized  equations  neglected 
by  Anderson  and  Rickayzen.  For  simplicity  we  treat  these  terms  only 
in  the 0  case .  In  the  equation  for  bt(t)  ,  the  terms  shown  in 
Pigs.  2(h)  and  2(i)  were  neglected.  They  contribute  the  factor 


k' 

Ifk'iYk'o^ 

^  ^  (i^.i) 

to  the  right  hand  side  of  (5. la)  in  the  limit 'q  -►  0,  while  the 
negative  of  this  factor  is  added  to  the  right-hand  side  of  (5*lb). 
The  exchange  scattering  vertex  shown  in  Fig.  5(a)  was  neglected 
in  the  equation  for  p  (^).  Its  contribution, 

Xu 

Z  I  - 

vanishes  as  q  -►  0  and  does  not  affect  the  energy  of  the  exciton 
states  in  this  limit.  The  inclusion  of  (4.l)  adds  the  term 


to  the  right-hand  side  of  (5- 6a)  and  the  negative  of  this  term  to 
the  rigtit-hand  side  of  (5«6b).  Introducing  the  variable 

4r' 

one  finds  the  M  0  exciton  states  satisfy  the  set  of  coupled 


equations : 
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■  —  ■  ■  /  j.  *7 


Rlm  ^  "y  Ak  _ 

^  f*  Lt.sif--‘ta^‘-  ‘ 

- 

(T  ^/<  (■t,iif--v-a„^ 


KchV,Z 


^lAk 


(4.5) 


Setting  the  determinant  of  the  coefficients  equal  to  zero,  one  finds 
the  dispersion  relation 

(k.6) 

or 

[ J—  AmiIL2L^  \(-l _ -1-  - 

i  x(j-%^)'^)\.Z~  f"  0-  V 

CAi^inxT  _  ^ 

/"  /-  ^ 


(4.7) 


for  the  energy  of  the  r_^  exciton.  The  modification  of  the  "5  0 

exciton  energy  given  by  (4-7)  is  shown  in  Fig-  6  for  g^  =  0.25  and 
is  seen  to  be  small.  A  new  type  of  excitation  follows  from  (4.7) 
for  g-  <  0,  that  is,  a  repulsive  rather  than  attractive  L-wave 
interaction  between  electrons.  The  energy  of  this  state  is  shown 
in  Fig.  8  as  a  function  of  -g,.  for  g  *  0.25*  From  the  form  of  the 

ii  u 

coherence  factors  entering  the  dispersion  relation  it  appears  the 
new  state  should  be  interpreted  as  a  bound  electron-hole  pair  in 
close  analogy  with  the  exciton  states  occurring  in  insulators.  This 
interpretation  is  consistent  with  the  fact  that  the  electron-hole 
interaction  is  attractive  when  the  corresponding  electron-electron 
interaction  is  repulsive.  Thus  the  electron-hole  exciton  arises 
solely  from  the  terms  neglected  in  the  Anderson-Rickayzen  equations. 

In  Sec.  V  we  consider  the  role  played  by  this  type  of 
exciton  in  the  absorption  of  electromagnetic  waves  of  frequencies 
lying  within  the  gap. 

V.  ELECTROMAGNETIC  ABSORPTION 

7/ 

Ginsberg,  Richards  and  Tinkham-^  have  measured  the  absorption 
of  infrared  radiation  in  bulk  seunples  and  the  transmission  through 
thin  films  of  several  superconductors  c  In  samples  of  lestd  and  tin 
they  have  found  a  precursor  absorption  existing  for  frequencies 
below  that  corresponding  to  the  energy  gap.  The  Rickayzen  form  of 
the  Anderson  equations  of  motion  has  been  applied  by  Tsuneto^^  to 
the  problem  of  the  surface  imi)edance  at  finite  frequency.  While  he 
finds  the  existence  of  a  precursor  absorption  due  to  the  creation 
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of  excitons,  his  results  predict  a  value  for  the  absorption  which 
is  an  order  of  magnitude  smaller  than  the  experimental  value.  The 


Tsimeto  analysis  does  not  include  the  corrections  to  the  Anderson- 
Rickayzen  equations  described  in  IV  which  give  rise  to  the  hole- 
particle  excitons.  In  this  section  we  calculate  the  infrared 
absorption  due  to  hole-particle  excitons  in  order  to  investigate 
how  this  absorption  may  modify  the  Tsuneto  results  and  to  see  if 
the  experimental  results  can  be  explained. 

In  order  to  calculate  the  absorption  we  must  extend  the 
equations  (5^1)  7^.  7,  and  7,  .  7,  to  treat  particle-hole 

excitons  with  a  finite  center  of  mass  momentum  yhq.  Once  again 
the  quasi -particle  conserving  operators  considered 

since  they  do  not  contribute  to  the  dispersion  relation  for  the 
excitons  nor  to  their  absorption. 

We  begin  by  defining  the  quantity  k* )  by 


With  this  definition,  the  non  quasi -particle  conserving  contribution 

to  the  equation  of  motion  (5* la)  for  7^  7 

K+qo  k^ 

of  the  terms  shown  in  Figs.  2(h)  and  2(1)  is 


due  to  the  inclusion 


The  contribution  to  equation  (5- lb)  is 

-c^(p-  V(^,K')+  VCk'^, 


(5.2b) 
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In  order  to  investigate  the  electromagnetic  absorption  due  to 
the  creation  of  excitons  we  consider  the  effect  on  the  equations  of 
motion  (3 -la)  and  (3.1b)  of  an  externally  applied  transverse  vector 
potential 

/4  .  (5.5) 

With  ot  =  (eH/2mc),  the  inclusion  of  the  vector  potential  and  the 
corrections  (5-2);  the  equations  (3.1a)  and  (3.1b)  are  rewritten  as 

-t  C„(p 

-I-  ^  L 

(5.^) 

and 

' i  h(l?.p^(p  -  4m  f)  -  4  ('p 

+  °<pC]lp)  tip)-Ul?-hp  exj>  [rc(uH>()tl. 


(SAb) 
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The  variables  A^(q),  B^(q),  and  have  been  defined  previously 

by  equations  (2.18),  It  is  the  three  coupled  integral  equations 
for  these  variables  which  are  the  basis  for  Tsuneto's  analysis 
of  the  electromagnetic  absorption.  With  the  inclusion  of  correction 
terms  due  to  the  presence  of  the  hole-particle  excitons,  the  two 
new  variables,  C  (J)  and  C  ('^),  of  equations  (5-2)  are  coupled  into 
the  equations  for  the  three  original  variables.  Tsuneto  has  shown 
that  the  cross  terms  involving  vanish  so  that  his 

absorption  is  determined  by  the  one  collective  variable  In 

our  analysis,  we  have  three  coupled  integral  equations  in  A^C^), 

C^(q)  and 

It  becomes  convenient  to  introduce  a  change  of  variables 
with  the  new  set  of  definitions i 

Tin  (p  =  Z  X 

'  Of  (5.5,) 

VM)x 

(5.5b) 


(5.5c) 
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We  also  note  that  D(^,lc+q)  +  D(lc-fqj'5)  - 

D('gM)  -  P(tt)  We  will  work  with  the  fol¬ 
lowing  three  coupled  integral  equations  in  ^k^^‘ 

\MCf  >{'  <p  4'f ^  ^  f  H ^ ^ 

(5.6a) 

^  ^  _/ _ 

IkCp-Z  X 

V 

(5.6h) 

Tp (p  =  ^ 

{t^iJWhppnp)  p^'.f)  ^(p-t' 

-■tMAy^',p/n(^',p/‘0?'pjlM'p}  pn. 
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We  again  use  the  decomposition  of  the  two  hody  potential  V(k*,'li) 
into  spherical  harmonics  as  described  by  equation  (5.7)-  This  allows 
us  to  define  the  three  quantities  A^(q),  '^2im^^^° 

4^a»i  J)  (fk^i'jo  j/k'i  ”/Il?  I 

(5.7a) 


(5.7c) 


Analogous  to  equation  (‘5.5c)  we  set 

’■!(«>  ■  W5>  •  (5-7d) 

With  these  definitions  euid  use  of  the  orthonormality  properties  of 
the  spherical  harmonics,  equations  (5-6)  may  be  reduced  to 
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Y*  Yc,  ulp 

YiM  -i)ij 

(5.8a) 


Xj^Cp  =  KX  X 

\^ti>MhpnC^',p ^  YX  Yen  Ajinp) 
-tcpp)>K(pp) (p ^ Y^ YxhZfi  (p 
~komCl?'p)7i^(pppi?'p)^  YJY  YmTah  (p 
-‘K^leimCP,pfCP,pd!(p)‘p  YYt  (fee,,  e-fn^-,)] J 


(5.3b) 
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(5^8c) 


We  must  calculate  the  correction  to  the  BCS  paramagnetic  current  due 
to  the  transverse  collective  excitations  which  are  included  in 
equations  (5«Q)«  Then,  following  the  method  of  Tsuneto,  we  determine 
the  surface  resistance  of  a  bulk  superconductor  for  frequencies  less 
than  the  gap. 


The  pareimagnetic  part  of  the  current  density  is 


(5.9) 
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By  the  methods  described  in  section  III,  we  take  matrix  elements 
of  equations  (5*4)  between  the  renormalized  ground  state  |  0^ 
and  the  state  |  l(q,a)^  containing  one  elementary  excitation. 

The  results  of  this,  together  with  our  e^qpressions  for 

+ 

allow  the  paramagnetic  current  density  to  be  rewritten  as 


-  -  Co{ 


\-Mt,pUp  kf 

ctp'Tp'lp 

-h  ^  A(p  ^(p  i*. 


(5.10) 

Since  the  odd  i  values  in  our  potential  expansion  (3.7)  can¬ 
not  couple  with  a  transverse  field,  we  consider  a  simplified  case 
in  which  the  only  terms  present  in  the  potential  are  and  V^. 
Fvtrthermore ,  the  presence  of  the  Kronecker  deltas  (due  to  the  trans¬ 
verse  nature  of  the  applied  field)  in  the  driving  terms  of  equations 
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(5.8)  indicate  that  we  need  only  consider  the  M  “  it  ^  cases  in 

solving  these  equations.  In  addition,  it  is  easily  shown  that 

+  + 

AgiCq)  =  -  and  T“^(q)  -  '^2,-1^^^* 

solve  equations  (5*8)  for  the  three  collective  coordinates 

and  T^^Cq).  We  reduce  our  notation  by  referring  to  these  as  A(q)  and 

T"(q). 

We  introduce  the  sums 


k'  (^r- 

i  ' 

-Li...  =  v^Z  , 


(5.11) 


(5.12) 


where  each  of  the  quantities  a,h,c,...  is  one  of  the  coherence  factors, 

y:  (^), /fioi  or  a(q)*k' .  With  this  notation,  the  three  equations  (5.8) 

+ 

for  A(q)  and  T  (q)  may  be  rewritten 


(5.13a) 


”  -h  ^n^tatrtnyT  =  fn^ 

(5.13b) 


/\  -hO-^  =  “W  'iy> 


(5.13c) 
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Vflien  the  angular  Integrations  are  considered,  the  orthonor¬ 
mality  properties  of  the  spherical  harmonics  reduce  equation  (5«13c) 
to 

Cl ^lhy^)^)'T  *0,  (5.14) 


Therefore^  a  consistent  solution  to  equations  (5*15)  may  obtained 
with  the  setting  of  T”  s  o,  since  the  quantity  (l  +  ^ 
non- vanishing  when  evaluated  at  points  of  the  dispersion  relation 
determined  by  A  and  t'*’  alone  o  The  two  equations  to  be  solved  are  now 


(5-15a) 


J  (5.15b) 

where  we  have  written  T**"  =  T  for  convenience  of  notation «  The  solutions 
are 

'  (5.l6a) 


(5.16b) 
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where 

?=  0-^ix^)(l4-S 

We  now  consider  the  calculation  of  the  surface  resistance.  In 
equation  (5.IO)  for  the  paramagnetic  current  density  the  first  term 
gives  the  BCS  paramagnetic  currents  The  collective  part  may  be 
written  as 

■joi  (fM  =  Z  -^^(p  ^(y^l  (% 

The  kernel  defined  by 

^Cd/  ^  Kcall  (p^)  ^ (p X  (5  18) 

and  can  be  easily  shown  to  be  given  by 


In  this  expression  the  quantity  ky  appears  as  an  approximation  to  (k| 
in  since  the  two-body  potential  is  assumed  to  vanish  except  for 
)k|  and  )k*|  near  ky. 
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To  evaluate  the  surface  resistance  we  assume  the  condition  of 
specular  reflection  of  the  quasi-particles  at  the  boundary  of  the 
superconductor.  We  should  not  expect  our  results  to  be  too  sensitive 
to  the  boundary  conditions.  For  example,  the  penetration  depth  com¬ 
puted  under  the  specular  reflection  condition  differs  little  from 
that  calculated  in  the  random  scattering  limit.  Also,  the  condition 
of  random  scattering  involves  mathematical  complexities  beyond  the 
scope  of  our  calculation. 

With  the  si)ecular  reflection  boundary  condition  the  instan- 

13/ 

taneous  field  inside  a  bulk  superconductor  is  given  by—' 


<^CTju)  =  _ i _  (5.20) 

2  2 

where  H(o)  is  the  magnetic  field  at  the  surface.  The  quantity  joT 
is  neglected  since  we  are  interested  in  the  infrared  region  where  the 
wavelength  of  the  electromagnetic  wave  is  large  compared  to  the 
penetration  depth  of  the  seunple.  For  our  rough  calculation,  we  also 
neglect  denominator,  so  that  K  21  K^,  the  ordinary  kernel. 

We  further  approximate  by  which  is  valid  for  ^^q  <  1.  The 

qiiantity  is  the  London  penetration  depth. 

The  rate  of  absorption  of  a  wave  of  frequency  o)  is 

LfCf,^)  •  t  (pa)]  di 

“  X,  L  Kca/ . 
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With  the  substitution  of  our  approximation  for  a(q,a))  into  (5-21), 
we  get 


Thus  we  are  interested  in  evaluating  the  imaginary  part  of 

K^oii(*^^^)  by  equation  (5*19)*  In  section  IIl(2)  e^qplicit 

expressions  could  not  be  obtained  for  the  integrals  involved  in 

A  (q)  and  r_^(q),  and  approximations  were  used  to  second  order  in 
In 

the  momentum  transfer  ^q  .  Similar  approximations  are  made  in  this 
section  for  the  evaluation  of  A,T  and  subsequently  Im 
With  these  approximations,  the  coupling  consteuits  defined  by  (5-8), 
and  X  =  (^/2A)  we  can  obtain 

Im  (W  X 

P)ic^O-kV  i- 


(5.23) 
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The  quantity  q(a))  is  the  solution  to  the  equation  P(q,a))  =  0  which 

gives  the  dispersion  relation  for  the  exciton  spectrum.  The  presence 

of  the  delta  function  term  indicates  the  absorption  by  the  system  of 

a  photon  of  frequency  oi  resulting  in  the  creation  of  an  exciton  of 

energy /fio)  and  center  of  mass  momentum  hq(a>).  Because  of  the  delta 

function,  the  integral  in  e:q?ression  (5*22)  for  R^Co))  is  easy  to 

o 

evaluate.  We  are  interested  in  obtaining  the  ratio-  of  R_(co)  to  the 

s 

resistance  of  a  normal  metal.  In  the  extreme  non-classical  (anomalous) 

limit,  and  with  the  boundary  condition  of  specular  reflection,  the 

ik/ 

surface  resistance  of  a  normal  metal  is  given  by — ' 


(u)=  {3)^(^^L/‘nOH§yA . 


(5.24) 


Prom  the  data  of  Ginsberg,  Richards  and  Tinkham  the  absorption 
in  the  gap  has  a  maximum  at  a  freqxiency  near  three-quarters  of  ?A. 

For  the  strong-coupling  metals  involved  in  their  experiments,  we  may 
take  a  value  of  g^  =  0.5*  In  Fig.  9  the  particle-hole  excitation 
energies  are  plotted  as  a  function  of  (-g^)  in  the  q  =  0  case  for 
this  value  of  g^.  We  see  that  em  excitation  given  by  x  =  ^  occurs 
for  a  value  of  g^  ~  -0.75*  With  these  values  of  g^  and  g^  the  ratio 
a>)  becomes 


6l 


X 


(5.25) 


This  ratio  has  heen  evaluated  as  a  function  of  x  =Jl(n/2A,  and 
the  results  are  indicated  in  Fig.  10.  We  see  that  our  results  do 
indicate  the  presence  of  absorption  due  to  the  creation  of  particle- 
hole  excitons.  We  do  not  observe  the  peak-like  structure  of  Tsuneto's 
results  since  our  approximations  are  not  valid  beyond  the  regions 
plotted.  However,  as  in  the  case  of  Tsuneto*s  work,  our  absorption 
is  at  least  an  order  of  magnitude  too  small  to  explain  the  observed 
experimental  data.  For  example,  with  equal  to  4  for  lead,  our 

results  indicate  an  absorption  ratio  of  about  0.005  at  x  =  0.825> 
while  the  observed  value  is  of  order  0.1.  Therefore,  we  can  conclude 
that  the  particle-hole  exciton  absorption  cannot  explain  the  lead  and 

7/ 

tin  precursor  absorption  observed  by  Ginsberg  et  al.,-'  within  the 
approximations  of  this  wesdt  coupling  calculation. 

It  should  be  noted  that  our  work,  as  well  as  that  of  Tsuneto, 
has  been  based  on  a  bulk  sample  calculation  while  the  experiments 
have  been  performed  on  thin  films.  The  wave  vectors,  q,  of  importance 
are  of  the  same  order  as  the  thickness  of  the  films.  Therefore,  an 
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improved  calculation  might  be  one  in  which  the  original  unperturbed 
wave  functions  are  not  three  dimensional  plane  waves,  but  ones  in 
which  the  finite  thickness  of  the  film  is  taken  into  account,  i.e., 
a  wave  function  which  is  an  infinite  plane  wave  in  two  dimensions 
and  sinusoidal  in  the  third  with  boundary  conditions  appropriate  to 
the  thickness  of  the  film  taken  into  account. 

VI.  CONCLUSIONS 

While  we  have  approximated  the  L  spherical  harmonic  of  the 

two-body  interaction  by  a  separable  potential v  V^(k,k*)  =  -  V.  for 

L  i-* 

1  ^kl  ^  l^k*l  zero,  otherwise^  in  general,  if  the 

potential  is  independent  of  crystallographic  orientation,  the  numbers 
L  and  M  remain  good  quantum  numbers  for  the  excitations  in  the  limit 
q  0.  For  a  nonseparable  potential,  ioC  . ,  if  V^(k,k*)  is  not  of 
the  form  ^yCk)  ),  there  may  be  more  than  one  exciton  state  for 

a  given  L  and  M.  While  the  excitons  should  give  a  negligible  contri¬ 
bution  to  the  specific  heat,  it  may  prove  possible  to  observe  the 
thermally-excited  odd  L  excitons  (spin  waves)  by  magnetic  resonance 
techniques.  The  precursor  infrared  absorption  observed  in  Pb  and  Hg 
cannot  be  explained  within  the  framework  of  a  bulk  seunple  calculation. 
It  would  be  interesting  to  carry  out  an  explicit  calculation  of  the 
absorption  coefficient  for  a  thin  film  geometry  in  an  attempt  to 
reconcile  the  difference  between  the  theoretical  predictions  and 


experiment. 
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Fig,  1  The  tvo  vertices  occurring  in  the  full  equations 
of  motion  for  bjlCq).  In  the  linearized  equations 
only  vertices  with  certain  values  of  p*  and  q* 


are  retained. 
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Fig.  1  The  two  vertices  occurring  in  the  full  equations 
of  motion  for  "b^Cq).  In  the  linearized  equations 
only  vertices  with  certain  values  of  p*  and  q^ 
are  retained. 
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Fig.  2  The  vertices  retained  in  the  full  linearized 

equation  of  motion  for  hMq).  Vertices  f^  h, 
and  i  were  neglected  by  Anderson  and  by  Rickayzen. 
The  particle-hole  excitons  are  obtained  only  if 


the  interactions  shown  in  h  and  i  are  included. 


6t 


Fig.  5  The  vertices  retained  in  the  full  linearized 
equation  of  motion  for  p,t(  q).  Vertices  a, 
e,  and  f  were  neglected  by  Anderson  and 
Rickayzen. 
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Fig.  k 


Fig. 


A  typical  diagram  retained  within  the  random- 
phase  approximation  to  P,  (q)  in  the  normal 
state . 


The  random-phase  approximation  to  the  screened 


interaction  line. 
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75 


Fig.  7  The  p-state  exciton  energy  as  a  function  of 

momentum  q  for  g^  =  0.25  and  g^  =  0.24  or  0.25* 

The  parameter  is  the  coherence  length 
*4 

~  10  cm.  Notice  that  the  exciton  states  are 
strongly  bound  only  for  q  ^  ^  ^q‘ 
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Fig.  9  The  energy  of  the  D-state  particle-hole  exciton, 
in  the  limit  of  q  0,  as  a  function  of  the 
D-wave  coupling  constant  g^  with  g^  =  0.50. 
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Fig.  10  The  ratio  of  the  surface  resistance  of  a 
superconductor  due  to  creation  of  D-state 
particle-hole  excitons  to  that  of  a  normal 
metal  in  the  extreme  anomalous  limit  cal¬ 
culated  for  the  D-wave  coupling  constant 


^2  =  ^0-75  and  g,  =  0.50. 
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Part  II 

Lifetime  Effects  in  Condensed  Fermion  Systems 
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I.  THE  ROI£  OF  LIFETIME  EFFECTS 
In  a  many  body  system  a  particle  in  an  excited  state  does  not 
remain  there  indefinitely  due  to  its  interacting  with  the  other 
elementary  excitations  of  the  system.  The  presence  of  this  lifetime 
can  be  taken  into  account  in  the  description  of  the  excited  state 
by  specifying  its  energy  as  E  =  E^  +  where  and  E^  are  both 

real.  From  this  we  see  that  a  well  defined  excitation  energy  only 
has  significance  when  the  state  is  long  lived;  i-e.,  when  E^ 

The  effects  of  finite  lifetimes  on  the  excitation  spectrum 
cannot  be  investigated  easily  using  an  equation  of  motion  approach 
or  a  variational  one  such  as  that  used  in  the  original—^  BCS  theory 
of  superconductivity  in  which  exact  eigenstates  are  obtained  for  a 
reduced  Heuniltonian  containing  no  damping  effects.  The  mathematics 
of  a  Green *s  function  formulation  provides  the  most  convenient  means 
for  investigating  the  lifetimes..  The  basis  for  an  analysis  of  this 
type  is  the  calculation  of  the  fermion  self-energy  diagram  of  Fig.  1 
whose  imaginary  part  is  related  to  the  damping  rate  The  solid  line 
represents  the  fermion  and  the  dotted  one  is  the  full  interaction 
line. 

2/ 

Migdal— '  has  evaluated  the  electron  self -energy  for  a  normal 
metal  with  the  interaction  line  representing  the  phonon  propagator. 
For  excitation  energies  (measured  relative  to  the  Fermi  surface) 

much  less  than  the  maximum  phonon  energy  JAay  he  finds  the  imeiginary 

o 

part  of  the  self -energy  to  be 


ail 


TrX^(?--K)  .  ‘^-1) 

In  this  e;q)res8ion,  is  a  dimensionless  parameter  of  the  order  of  l/2 
related  to  the  coupling  strength  between  the  electrons  and  phonons. 

For  energies  €.  >>tta)  ,  Migdal  finds  that  the  deunping  rate  due  to  real 

K  0 

phonon  emission  saturates  out  to  a  constant  value »  For  the  range  of 
excitation  energies  larger  than  about  ten  times  the  maximum  phonon 
energy,  damping  due  to  hole-electron  pair  production  becomes  important 
and  has  been  obtained  by  Quinn  and  Ferrell-^  by  allowing  the  inter¬ 
action  line  in  the  self -energy  diagreun  to  represent  the  screened  Coulomb 

interaction.  For  values  of  e,  <  E  ,  the  Fermi  energy,  they  obtain 

K  r 


2  1  /2 

where  -  (Uicne  /m)  *  is  the  plasmon  frequency. 

In  the  case  of  a  transition  to  a  superfluid  state,  the  effect 
of  damping  on  the  required  energy  gap  and  the  transition  temperature 
is  closely  related  to  the  relative  roles  the  phonon  and  Coulomb  inter¬ 
actions  play  in  constructing  a  criterion  for  the  possibility  of  a 
transition. 

In  the  original  BCS  theory  the  general  criterion  for  super¬ 
conductivity  is  that  the  attractive  electron-phonon  interaction  as 

4  / 

calculated  by  Bardeen  and  Pines—  dominate  the  Coulomb  interaction  for 


those  matrix  elements  which  are  important  in  the  superconducting  wave 


function.  The  specific  model  treated  ty  BCS  is  one  in  which  the 
two-body  interaction  is  taken  to  be  the  separable  potential 


-  Vkk'  = 

-p>r 

—  -  0  ,  otherwise. 

In  this  expression,  and  ai-'e  tbe  Coulomb  and  phonon  interactions 
and 

AY\r\ 

5/ 

The  BCS  criterion  has  been  further  investigated  by  Pines—  who 
considers  the  sign  of  the  interaction  (V  +  V  )  =:  -V  ,  in  the 

C  pn  KK  KK. 

zero  excitation  energy"  limits  l.e.,  when  both  |1?  land  lie  *1  are  taken 
equal  to  the  Fermi  momentum  k^.  In  this  model  the  ionic  charge  is 
considered  to  be  a  continuous  fluid  whose  vibrational  spectrum  pro¬ 
vides  the  frequency  distribution  of  the  phonons  used  in  the  problem. 
With  this  model.  Pines  finds  V  ,  will  be  positive  at  the  Fermi  surface 

iUv 

if  Umklapp  processes  are  considered  (which  are  calculated  without 
including  periodicity  effects)  Using  this  criterlony  that  the  net 
interaction  must  be  positive  at  the  Fermi  surface;,  predictions  are 
made  as  to  which  elements  in  the  i>eriodic  table  ought  to  be  sui>er-’ 
conducting.  In  the  case  of  those  elements  for  which  the  neglect  of 
periodicity  in  calculating  the  matrix  elements  may  be  valid,  certain 
qualitative  results  are  obtained  which  agree  with  experiment. 
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Another  investigation  into  the  criterion  for  superconductivity 
has  heen  made  by  Bogoiiubov^,^  and  Abrikosov  and  Khalatnikov.— ^  In 
their  treatments,  the  Coulomb  and  phonon  interactions  are  considered 
explicitly  in  order  to  determine  the  relative  roles  played  by  these 
two  interactions  in  the  criteriono  In  their  method^  the  attractive 
potential  due  to  the  phonon  interaction  only  is  taken  into  account  by 
inclusion  of  a  separable  potential  (l.5)>  and  is  characterized  by  the 
positive  coupling  constant  g  -  N(o)V  >  0.  The  q^uantity  N(C/  is  the 
density  of  Bloch  electron  states  of  one  spin  at  the  Fermi  surface c  The 
Coulomb  interaction  is  accounted  for  by  a  constant  repulsive  potential 
acting  within  the  momentum  band,  which  is  larger  than  the  range, 

over  which  the  attractive  potential  acts.  This  corresponds  to 
the  fact  that  the  Coulomb  interaction  is  screened  at  distances  of  the 
order  of  the  lattice  spacing  and  will  be  explained  in  greater  detail 
below.  The  repulsive  potential's  coupling  constant  is  denor^ed  by 
g*  -  -N(o)V*  >  0.  With  these  definitions  the  gap,  A,  is  the  solution 
of 


/ 


— * 


(1.4) 


Prom  this  equation  we  see  that  the  criterion  for  super? or ductivity  in 
this  model  is 


; 


> 


(1=5) 


replacing  the  condition  g  >  g'  given  hy  the  separable  potential  model 
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of  BCS.  The  role  of  the  Coulomb  interaction  is  thus  reduced  by  the 

factor  ln(uj  */“  )  which  is  typically  about  5  if  hco’  •:l  E_..  Furthermore, 
c  c  c  r 

when 


(1.6) 


the  net  potential  may  be  everywhere  repulsive  in  mcmentum  space  and  still 
lead  to  a  superconducting  state.  This  result  may  be  more  clearly  under¬ 
stood  if  we  consider  the  relationship  between  an  interaction  in  coordinate 
space  and  its  momentum  space  Fourier  transfcrra.  For  the  limiting  case 
of  a  constant  potential  having  infinite  extent  in  momentum  space,  the' 
coordinate  space  interaction  takes  the  form  of  an  infinitely  sharp  delta 
f‘anction.  As  the  e^rtent  of  the  interaction  in  momentum  space  is  shortened, 
the  corresponding  delta  function  broadens  increasing  the  interaction's 
extent  in  coordinate  space.  We  can  now  apply  this  notion  to  the 
Bogoliubov  potential  in  which  the  Coulomb  Interaction  has  a  greater 
extent  in  momentum  space  than  the  phonon  interaction.  With  ^ 
the  coordinate  space  extent  of  the  screened  Coulomb  interaction  (the 
order  of  a  lattice  parameter)  is  less  than  that  of  the  phonon  inter¬ 
action.  From  the  BCS  theory  we  know  nhat  the  superconducting  state 
is  characterized  by  the  formation  of  behind  pairs  of  electrons  whose 

relative  coordinate  wave  function  has  er  ext.en+  of  10  cm-.  Since 

-8 

a  lattice  parameter  is  roughly  10  cm^  the  bound  pair  of  electrons 
rarely  experiences  the  Coulomb  potential^  although  its  effects  may  not 
be  totally  neglected  since  it  is  of  large  magnitude.  The  bound  state 
is  mainly  determined  by  the  attractive  phonon  interaction  even  though 
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its  average  magnitude  may  be  less  than  that  due  to  Coulomb  effects. 

In  this  way  a  net  repulsive  interaction  in  momentum  space  can  lead  to 
a  superconducting  transition.  Of  course,  when  the  range  of  the  two 


interactions  is  the  same  so  that  cn^  *  -  the  Bogoliubov  criterion 


reduces  to  that  of  Pines  and  the  model  solved  by  BCS.  One  difficulty 
in  greatly  reducing  the  Coulomb  effects  is  that  the  critical  tempera¬ 
ture  increases  so  that  most  metals  might  be  e2q>ected  to  exhibit  super¬ 
conductivity. 

O  j 

Bardeen—'^  has  suggested  that  lifetime  effects  due  +-o  the  Coulomb 
excitation  of  particles  out  of  the  Fermi  sea  could  be  used  to  determine 
the  cutoff , *,  for  the  Coulomb  interaction.  When  the  energy  ^  is 
sufficiently  large,  the  excitation  may  decay  so  rapidly  that  it  is 
not  well  defined.  For  a  rough  estimate  of  Coulomb  damping  effects  we 
can  use  the  Quinn  and  Ferrell  result  of  equation  (ln2).  From  this 
equation  we  see  that  for  ^  Ep,  the  imaginary  part  of  the  self  energy 
is 

j  Cf:  .  ^ 


As  we  will  see  in  the  next  section,  damping  begins  to  affect  +he  gap 

equation  when  E^(k)  and  are  equal.  Although  equation  (21)  does  not 

hold  for  ^  result  (1.7)  indicates  that  the  desired  equality 

would  not  tedce  place  until  c,  >  E  .  But  for  these  values  of  c,  the 

K  r  K 

matrix  elements  of  the  screened  Coulomb  interaction  are  small,  indi¬ 
cating  that  Coulomb  danqping  effects  play  a  small  role  in  determining 

the  cutoff /fto)  *. 

c 
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This  result  shows  that  the  inclusion  of  Coulomb  deunping  effects 
cannot  reduce  the  logarithm  term  in  the  Bogoliubov  result  (1.5)'  This 
reduction  is  necessary  if  the  role  the  Coulomb  interaction  plays  in 
establishing  a  criterion  is  to  be  increased.  On  the  other  hand^  as  a 
consequence  of  (1.5)^  most  metals  would  be  expected  to  be  superconducting. 
We  can  only  conclude  that  a  careful  examination  of  periodicity  effects 
in  the  electron-phonon  interaction  would  be  necessary  in  order  to  estab¬ 
lish  a  satisfactory  criterion  for  superconductivity^  This  problem  is 
not  considered  in  this  thesis. 

Below,  we  outline  a  formalism  for  working  with  damping  effects  in 
a  fermion  system  and  consider  the  consequences  of  certain  assumed  forms 
for  this  deunpingo  We  also  consider  the  problem  of  He^  for  which  no 
superconducting  transition  has  been  observed  in  experiments  down  to 
T  >  5  X  lO’^^K,  although  theoretical  predictions  neglecting  damping 
have  predicted  a  critical  temperature  an  order  of  magnitude  larger  than 
this  value.  We  find  that  the  damping  of  excited  states  will  tend  to 
reduce  the  net  attractive  interaction  necessitating  a  corresponding 
lowering  of  the  critical  temperature  > 

II.  EVALUATION  OF  THE  SELF -ENERGY  DIAGRAM 

In  calculating  the  self  energy  diagram  of  Fig-  1  ve  will  follow 

9  ^  10  ’ 

most  closely  the  formulations  of  Nambu-  and  Eliashberg  —  It  is 
useful  to  introduce  the  two-component  form  for  the  electron  wave 
fields 


(2.1) 
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along  with  the  Pauli  spin  matrices 


i;. 


(2.2) 


The  operators  C  and  C  are  the  usual  electron  creation  and  annihila- 
^  pa  po 


tion  operators  whose  only  non-vanishing  anti -commutators  are 


1  C p<r^  Cp‘(r‘  1  STcnr'  ^ 


(2.5a) 


from  which  we  see  that 


(2-5b) 


In  evaluation  of  Feynman  diagrams  within  the  new  notation  of  ^2.l) 
and  (2.2),  all  the  usxoal  rules  hold  with  the  ^*s  taking  th#e  roles  of 
the  c*s  of  the  old  notation.  The  only  new  rule  is  the  addition  of  a 
factor  of  Tj  at  each  vertex<  This  can  be  seen^  for  example  by  con¬ 
sidering  the  feunilar  kinetic  energy  and  electron -phonon  interaction 
Hamiltonians 


and 


Hk-  -  Cpo-  Cp<r 

f,o-  ^ 


1 2 .  iia ) 


(2.1rt)) 


r2_2 


where  €  =  -  E_,  g  is  the  strength  of  the  interaction,  and  a 

p  ^  r  q  q 

is  a  bare  phonon  operator.  In  our  new  notation  equations  (2«4)  become 


Hic=Z.^p[  5f^^,iFp+0, 


(2.5a) 
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(2.5b) 


If  we  consider  the  Fermi  sea  as  our  vaeuiam^  1  0^  .  we  define 
the  Green's  function  for  free  electrons  (H  -  hy 


J  }  (2.6a) 


where  T  is  the  time  ordering  operator.  Taking  the  Fourier  transform 
of  (2.6a),  it  follows  that 


(2.6b) 


where  p  -  (?»p  )•  The  Green's  function.  G  ,  for  the  total  Hamil‘*’.onian 
H  *  +  H*  is  related  to  the  self  energy,  I  and  the  lowest  order 

IV  p 

Green's  function  (2.6b)  through  the  Dyson  equation 


(2.7) 


The  most  general  form  for  would  be  one  having  components  along  all 
four  orthogonal  axes;  i*e.. 


Zp=-t(^4r  + 


(2.8) 


However,  the  first  order  Hamiltonian  (2.5a)  is  invariant  under  rota¬ 
tions  in  a  plane  orthogonal  to  the  unit  matrix  and  7,.  Therefore,  we 

5 

con  set  0  in  (2.8)  and  take 


(2.9) 
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This  equation,  together  with  (2.6b)  and  (2^7),  gives 

po  ^  (Xp-^Xf,)  tj  *C,  ^ 

where 

4=  /  -  c^/pj. 


(2.10) 


(2.11) 


The  Dyson  equation  connecting  the  self-energy,  E  ^  with  the 

P 

exact  Green’s  function  for  the  pro’uiem  and  the  total  vertex,  F  ..is 

pp ' 


(2.12a) 


where  ^ppi  the  propagator  for  the  interaction-  The  self-energy  is 
generally  obtained  from  a  perturbation  expansion «  In  the  case  of  an 

Q 

electron-phonon  interaction,  ^ppi  phonon  propagator  and  ip 

represents  the  coupling  constant.  For  this  case  Migdal—'^  has  snown 

that  r  ,  Cl  to  order  (m/M)  where  M  is  the  ion  mass.  If  the 

P  *P  P  'F 

vertex  function  for  einy  interaction  is  approximated  by  its  zeroth  order 
value,  this  is  equivalent  to  summing  a  perturbation  series  for  the  self* 
energy  only  Involving  diagrams  in  which  no  two  interecticn  lines  cross. 
The  result  of  this  eunounts  to  evalizating  the  self -energy  diagram  of 
Fig*  1  in  which  the  solid  line  represents  the  full  Greenes  function  for 
the  problem  and  the  dotted  one  the  full  interaction.  For  a  general 
potential,  ^ppi>  approximation  to  equation  (2- 12a)  Is  expressed  by 


(2.12b} 
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where  the  T*s  appear  for  the  reasons  discussed  previously.  Substituting 
(2*i0)  into  (2»12b)  and  considering  the  component  of  the  result,  we 
obtain  the  gap  equation 


_ (2.15) 

The  integration  over  p*  can  be  carried  out  by  the  considerations  of 


Eliashbergo 


The  Green’s  fimction  G  ..  when  considered  as  a  fimction 

P' 


of  Pq,  is  analytic  in  the  upper  (lower)  half  plane  for  >  (<)  Oc 
Therefore,  the  i^pi  berm  may  be  dropped  in  the  integrand  when  the 
integration  is  taken  along  the  contour  of  Fig,  2.  We  will  now 
assume  that  the  potential  is  real  and  independent  of  and  p^. 

With  this  condition,  consideration  of  the  complex  function  Gfp.z) 
which  coincides  with  G(p,Pq  )  on  the  upper  (lower)  half  plane  for 
>  (<)  0  allows  the  contour  to  be  deformed  to  of  Figo  5 
Since  the  values  of  iG(p,z)  on  the  opposite  sides  of  the  cut  (O,®) 
are  complex  conjugates,  we  obtain 


oo 

'  J  iff  V  Sr) 


(2-1^) 


where 


£,■  = 


(2.15) 


Since  the  integreuid  is  symmetric  in  p^'  the  lower  limit  in  (2.ll4)  may 
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be  replaced  by  -oo  and  the  result  multiplied  by  1/2.  Assuming  that 
the  analytic  continuation  of  Into  the  lower  half  plane  for  p^  >  0 
has  a  simple  pole  corresponding  to  an  elementary  excitation  of  energy 


(2.16) 


tue  integral  eqiiation  becomes 


(2.17) 


In  this  expression  the  quantity  related  to  through  (2cl2)^ 

satisfies  an  integral  equation  similar  to  (2.15)  obtained  by  considering 

the  unit  matrix  element  component  of  equation  (2. 8)  and  gives  the 

normal  state  self-energy  in  the  limit  of  A  0»  The  imaginary  part 

P 

of  gives  rise  to  damping  effects  in  the  energy  gap  equation.  When 
^  is  set  eqxial  to  zero  (Z^  -  1;,  the  integrand  of  (2  >17)  becomes  real 
and  the  equation  reduces  to  that  of  BCS. 


Ill,  DAMPING  EFFECTS  FOR  ENERGY  INDEPENDENT  V  , 

P'P 

We  will  investigate  the  effects  of  damping  on  equation  (2.17)  by 
considering  various  forms  of  the  function  The  real  part  of  ^ 

will  not  be  needed  since  it  serves  to  renormalize  the  excitation  ener¬ 


gies  Cp  and  the  interaction  ^ppt  imaginary  part  of  ^ 

which  provides  for  the  damping  effects.  We  consider  a  general  form 
for  in  which  the  dsuoaping  is  proportional  to  a  power  of  p^i  i.e., 


-If-  t  p<  IppI 


(5.1) 


95 


It  should  be  noted  here  that  our  assigned  form  for  ^*8  imaginary 
part  is  not  entirely  satisfactory  since  this  quantity  must  vanish 
for  |p^|  Therefore,  for  n  «  0  in  (5*1)  we  may  onljj  consider, 

the  0  situation.  For  u  »  1  the  approximation  is  much  better, 

and  for  the  case  where  n  ~  2,  the  error  introduced  is  entirely 
negligible.  For  n  ^  0,  substitution  of  (5-1 )  into  (2.1?)  gives 


where  is  the  solution  to  -  E  « 

^0  ^0  p* 

We  will  consider  the  three  cases  n  -  0,  1,  and  2. 


(1)  n  -  0 


The  n  =  0  case  corresponds  to  constant  damping.  For  the 
separable  two  body  interaction 


(5.5) 


Vpp'  - 

setting  Ap  =  0  in  the  n  =  0  integral  equation  allows  us  to  determine 

the  minimum  value  of  a  such  that  no  transition  to  the  superconducting 

state  occurs;  i.e.,  =  0°K,  Although  equation  (5.^}  does  not  hold 

for  n  0,  it  is  easy  to  show  that  the  A  0  equation  for  this  case 

P 

is 


/ 

mow 


) 


(5.4) 


where  N(0)  is  the  density  of  Bloch  states  of  one  spin  per  unit  energy 


at  the  Fermi  surface.  The  solution  to  (3«^)  gives  the  critical 
value  for  a, 

<SX^[-lMo)V]  ,  (5.5) 


This  result  agrees  with  that  of  Svihl 


11/ 


(2)  n  -  1 

With  the  separable  potential  of  case  (l)^  the  gap  equation  (5 <>2) 


gives 

J-  ^  f  di.  ! 

A^(P)V  ^ 

The  solution  is  of  the  BCS  type  with 


(5.6) 


Ap=  ^COq,  l£pl  —'hlAi  ' 

A  =  0  ,  othervir:.e. 

p 

Thus,  the  inclusion  of  linear  damping  modifies  the  BCS  solution  only- 
through  the  introduction  of  a  reduced  effective  coupling  constant 
N(0)V/(1  +  a^). 


(5)  n  -  2 

For  the  case  of  quadratic  damping,  the  T  •=-  O^K  gap  equation 


is 


From  the  result  of  (l«7)  ve  have  seen  that  the  screened  Coulomb 
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interaction  leads  to  l/a  >  E_  so  that  the  effects  of  Coulomb  damping 

c 

on  equation  (3*8)  are  negligible  since  the  screened  Coulomb  matrix 

elements,  V  ,jare  small  for  energies  €  >E  .  Therefore,  inclusion 
PP  p  F 

of  Coulomb  damping  cannot  serve  to  decrease  the  extent  of  the  Coulomb 
interaction  in  (5*8)o  Consequently,  the  parameter  in  the 
Bogoliubov  cri'cerion  (1.5)  for  a  superconducting  state  in  metals  is 
essentially  unchanged  from  its  value  computed  in  the  absence  of 
damping. 

For  He^,  damping  effects  are  important  in  determining  and 
the  transition  temperature,  T^,  for  a  possible  superfluid  state.  Due 
to  the  hard-core  potential,  the  integral  equation  must  be  solved  in 
coordinate  space.  However,  we  can  obtain  axi  estimate  for  the  reduction 
of  due  to  damping  effects  by  considering  an  effective  separable 
potential  given  by  (3»3)»  With  this  potential,  the  transition  tem¬ 
perature  is  given  by 

/^c. 

J  _  (5  „) 

-Jo  _/J 

where  6  -  (kT 

We  will  now  calcxilate  an  approximate  value  for  a  For  the 
thermal  conductivity  of  He^,  Abrikosov  and  Khalatnikov^'^  obtain  the 
theoretical  expression 

/C  =  -f-  I  r 

7”  L  ^  (^f) 

where  p^  is  the  Fermi  momentum  emd  m*  is  the  effective  mass  of  a  He^ 


j  (5.10) 
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atom.  The  quantity  a)(0,cp)  represents  the  transition  rate  for  the 
scattering  of  a  quasi-particle  and  is  related  to  the  scattering 
eunplitude,  f(0,<p),  "by  the  usual  "golden  rule"  formula 


(5.11) 


Since  we  have  chosen  an  angularly  independent  separable  potential  it 
is  consistent  to  consiuer  as  constant  in  the  expressions  (5*10) 

and  (5*ll)*  With  this  assumption,  performing  the  angular  average  in 
(5*10)  gives 


C(r= 


(5.12) 


At  this  point  we  consider  the  above  quoted  work  of  Quinn  and 
Ferrell,^  in  which  the  imaginary  part,  excitation 

energy  is  obtained  for  electrons  interacting  under  the  screened 
Coulomb  potential.  Their  results  may  be  used  for  the  He^  problem 
if  the  screened  Coulomb  matrix  elements  are  replaced  by  the  scattering 
eunplitude,  f.  This  straiglitfcrward  generalization  gives 


3  2 

where  n  is  the  He  particle  density.  Substituting  f  obtained  from 

(5.11)  and  (5-12)  into  (5.15^  g'-ves 


(5.14) 


Use  Is  now  made  of  the  e3q>eriinental  work  of  Anderson,  Salinger  amd 
Wheat  ley^^  Mho  obtain 
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(3.15) 


With  this  value  for  the  thermal  conductivity^  equation  (5»1^)  becomes 

E,  cf.)  =  i/--  ^ 

22  /  3 

The  values  m*  =  2.82m  and  n  -  1.64  x  10  particles/cm^  give  the 
deunping  constant 

a  “  4.27  X  lO^^lergs)  ^  (Jol?) 


5  0-1 

When  compared  to  the  He  Fermi  temperature  of  about  5  K,  a  “  lo77  > 

indicating  that  the  imaginary  and  real  parts  of  the  excitation  energy 

are  equal  at  e  ~  E-/2.  Having  arrived  at  a  value  of  a>  ve  may  go 
P  ^ 

on  to  compute  T^. 

The  integral  of  equation  (3o9^  is  approximated  by 


/V(o)^ 


4 


-f 


i. 

/ _ 

-i] 


(3.18) 


For  this  approximation  to  be  reasonable  we  need  a  \alue  of  for 
2  2 

which  ^  r  «  1  for  €  <  In  this  way  damping  effects  are 

negligible  in  the  c  range  of  I.  and  hyperbolic  tangent  may  be 
replaced  by  unity  in  We  note  that  for  =  l/lO  a  the  condition 

of  is  satisfied.  Furthermore,  if  we  take  T^  -  0.05^K  for  the 
critical  temperature  without  damping  as  calculated  by  Emery  and 
Sessler,— our  choice  of  ^  ^  gives  tanh(p^  ^,^2)  0.9I«  Since 
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our  result  for  T  will  be  less  thsua  T°,  setting  =  l/lO  a  also 

C  C  X 

satisfies  the  requirement  of 

The  integral  is  identical  to  that  arising  in  the  BOS 
critical  temperatiore  equation.  Their  numerical  evaluation  gives 


(5.19) 


For  =  1/10  a,  the  upper  limit  in  is  essentially  infinite^  so 
that  we  get 


(5.20) 


With  these  approximate  values  for  and  in  equation  (5.I8),  the 
critical  temperature  with  inclusion  of  damping  effects  is  given  by 


=  (^soyOJi)  i,  e/p  [-  vj , 


(3.21) 


The  ratio  of  T  to  T°  is  then 
c  c 

=  C^s'o)^^  (£,  / 'IqLOd), 


(5.22) 


For  €  ^  =  1/10  a  =  Ep/20  and  ~  E^, 

-  0,32^  (3-25) 

indicating  that  the  effects  of  damping  in  our  simplified  treatment 
would  reduce  the  critical  temperature  by  a  factor  of  about  three. 

For  the  Emery  euid  Sessler  value  of  T°  =  0.05°K,  damping  effects 
reduce  the  critical  temperature  to  T^  0.0l6°K.  Ejqieriments  have 
been  performed  at  temperatures  T  >  5  x  10"^°K  without  the  appearauice 
of  a  superfluid  transition,  which  might  indicate  that  damping  is  not 
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as  important  in  determining  as  we  had  hoped.  However,  the  reduction 

of  T°  hy  a  factor  of  three  could  be  important  if  a  lower  valtoe  of  T° 
c  c 

were  ol)tained.  We  should  also  add  at  this  point  that  our  work  was 

based  on  a  simplified  model  in  which  the  potential  was  assumed  to  be 

of  the  S-wave  type;  i,e.,  V  ,  contains  no  angular  dependence.  On  the 

other  hand,  the  work  of  Emery  and  Sessler  is  based  on  the  assumption 

of  a  D-wave  term  in  the  potential.  Therefore,  it  might  not  appear  valid 

for  us  to  use  their  result  of  Ci  Oo05°K  in  our  calculationo  We  have 

c 

seen,  however,  that  the  effect  of  damping  is  essentially  to  reduce  the 
amount  of  phase  space  over  which  the  potential  acts.  For  this  reason 
the  ratio  would  not  be  expected  to  vary  significantly  with  the 

choice  of  angular  dependence  for  the  potential,  although  in  a  more 
exact  calculation  of  a  D-wave  term  would  explicitly  be  included 
in  the  potential  in  the  energy  gap  equation. 
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Fig.  1 


Fig.  2 


Fig.  3 


Self-energy  diagreim  used  to  evaluate  the  integral 
equation  for  the  energy  gap 


Contour  C^  along  which  the  integration  is  carried 
out  in  the  energy  gap  equation  (2.13)- 


For  an  energy  independent  potential,  V  ^  ,  the 
contour  0^  may  be  deformed  to  leading  to  the 
form  (2.1^)  for  the  energy  gap  equation. 
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